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Chandler wobble: Stochastic and deterministic 
dynamics 


Alejandro Jenkins 


Abstract We propose a model of the Earth’s torqueless precession, the “Chandler 
wobble”, as a self-oscillation driven by positive feedback between the wobble and 
the centrifugal deformation of the portion of the Earth’s mass contained in circu¬ 
lating fluids. The wobble may thus run like a heat engine, extracting energy from 
heat-powered geophysical circulations whose natural periods would otherwise be 
unrelated to the wobble’s observed period of about fourteen months. This can ex¬ 
plain, more plausibly than previous models based on stochastic perturbations or 
forced resonance, how the wobble is maintained against viscous dissipation. The 
self-oscillation is a deterministic process, but stochastic variations in the magnitude 
and distribution of the circulations may turn off the positive feedback (a Hopf bifur¬ 
cation), accounting for the occasional extinctions, followed by random phase jumps, 
seen in the data. This model may have implications for broader questions about the 
relation between stochastic and deterministic dynamics in complex systems, and the 
statistical analysis thereof. 


1 Introduction 

According to Euler’s theory of the free symmetric top, the axis of rotation precesses 
about the axis of symmetry, unless both happen to align IT]- Applied to the Earth, 
this implies that the instantaneous North Pole should describe a circle about the 
axis of symmetry (perpendicular to the Earth’s equatorial bulge), causing a periodic 
change in the latitude of any fixed geographic location (an effect known as “polar 
motion” or “variation of latitude”). The angular frequency of this free precession is 
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h-I 

(Oeu = —J —^23 , (1) 

where I = I\ = h and are the Earth’s moments of inertia, and Qj, is the compo¬ 
nent of its angular velocity along the axis of symmetry, which is also the principal 
axis with moment of inertia I^. Many textbook treatments of rigid body dynamics 
cover this in detail; see, e.g., . (Some authors refer to the corresponding terrestrial 

motion as the Earth’s “free nutation”; see El.) 

The value of I/{I^— /) can be deduced from the periods and amplitudes of the 
Earth’s forced precession and nutation, driven by tidal torques El- The Eulerian 
period of the Earth’s free precession then comes out to InjcOEu = 306 days (about 
10 months), but Chandler’s careful observations, first published in 1891, established 
that the actual period is about 14 months EUSlIbl. This motion is therefore known 
as the “Chandler wobble”. A modern estimate Q of its angular frequency is 

_ 2n 

“ 433.0±1.1 (la) days ' 

In 1892, Newcomb explained the discrepancy between (Uch and oJeu as a conse¬ 
quence of the Earth not being quite rigid. ID 

Newcomb also opened a debate, which has raged since, on how the Chandler 
wobble is maintained against the internal viscosities that would damp it away. Pro¬ 
posed solutions to this puzzle have, for the most part, relied either on stochastic 
perturbations to the Earth’s mass distribution Esni, or on marine and atmospheric 
circulations that might force the wobble near its resonant tOch. muiminiiTi 
Here we propose a new model of the Chandler wobble as a weakly non-linear 
self-oscillation. Unlike a resonator, a self-oscillator maintains periodic motion at 
the expense of a power source with no corresponding periodicity ifTSl fl^ . A self¬ 
oscillation modulates the driving force acting on it, establishing positive feedback 
between it and the external power source. On the many names and guises of self¬ 
oscillators, see ED. One striking instance of such behavior was provided by the 
swaying of London’s Millennium Bridge when it opened in 2000, driven by the 
synchronized lateral motion of the pedestrians in response to the bridge’s own os¬ 
cillation. E3 

In our model, the Chandler wobble is powered by geophysical fluid circulations 
whose natural periods are unrelated to Och. Due to the centrifugal force of the 
Earth’s rotation, the wobble itself can modulate those circulations. That a dynamical 
delay in the adjustment of a surrounding flow to a solid’s displacement can excite 
and maintain the solid’s vibration is an idea dating back to Airy’s conceptualiza¬ 
tion of the action of the vocal cords ED, based on Willis’s pioneering research on 
the mechanism of the larynx ED. In the case of the Chandler wobble, variations 
of the magnitude and geographical distribution of the circulations may turn off the 
positive feedback. This suggests an explanation of the wobble’s extinctions, which 
occur rarely and are followed by re-excitation with a random phase jump, without 
obvious connection to major geophysical events. 1^ 
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2 Precession and deformation 

Newcomb’s simplified treatment in IS) of the Earth’s deformability and its connec¬ 
tion to the observed coch lends itself to a physically intuitive formulation of our 
self-oscillatory model. We therefore begin by summarizing Newcomb’s argument. 
For a full, rigorous treatment of this subject, see 0. 

Let £l be the Earth’s angular velocity and Jcs be the unit vector along the axis of 
symmetry, with moment of inertia is I^. Let I be the value of the moments of inertia 
along any two axes orthogonal to each other and to jc 3 . Since the Earth is fiattened 
at the poles, > I. As shown in Fig. [1] for a rigid Earth the instantaneous North 
Pole R (the intersection of Q with the Earth’s surface) would describe a circular 
trajectory (polhode) around the symmetry pole S (the intersection of Xt, with the 
Earth’s surface), with angular velocity Weu given by Eq. O- 

The centrifugal force generated by the planet’s rotation implies that a displace¬ 
ment of R tends to deform the planet by shifting its equatorial bulge, causing S to 
move towards R. If the Earth were simply fiuid, this adjustment would be complete 
and nearly instantaneous, so that S would always coincide with R and there would 
be no wobble. For an elastic Earth, the adjustment is partial and the planet may 
wobble, albeit more slowly than in the rigid case. 

Let O be the center of the polhode for the elastic Earth (i.e., the average location 
of R over a complete period of the wobble) and let us work in a geographic frame of 
reference, centered at O. As shown in Fig.|2] S describes its own circular trajectory, 
with fixed 

kw = OS/OR (3) 

related to the speed of rotation £1 and to the planet’s elastic modulus. The parameter 
kw lies between between 0 (rigid case) and 1 (fiuid case). 


Fig. 1 Free precession of 
a rigid Earth. The principal 
axes jci, 2,3 define a non- 
inertial body frame, being 
the Eaifh’s axis of symmetry. 
These principal axes move 
with respect to the inertial 
space frame (for which only 
the axis Z is shown). S is 
the symmetry pole and R the 
instantaneous North Pole. In 
the body frame, £i precesses 
regularly along the circle SR, 
with angular velocity £Oeu . 
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The centrifugal deformation preserves the Earth’s ellipticity, leaving the magni¬ 
tudes of I and It, unaffected. Newcomb’s insight was that the planet’s internal forces 
of cohesion will therefore make R = (x,y) precess instantaneously about S = 
with the same Weu of Eq. O: 

rx=-a)Eu(y-ys) 

\y = a)Eu(x-x5) 

By Eq. (|3]l, X 5 = kwX and ys = kwj, so that Eq. (HJi reduces to 

^Cny ^ ( 5 ) 


with 

®Ch = ®Eu (1 ) 

the angular velocity of the precession of R around O. The parameter 

®Eu 433 


( 6 ) 

(7) 


measures the Earth’s deformability and is consistent with the Love number kj com¬ 
puted from the magnitudes of terrestrial and oceanic tides. Euim 

This simple analysis reveals that the dynamics of a deformable spinning body is 
subject to a feedback: the displacement of R (the wobble) affects the displacement 
of S (due to centrifugal deformation), while the displacement of S (the deforma¬ 
tion) in turn affects the displacement of R (due to precession). If some of the mass 
of the planet is in the form of circulating fluids (with their own kinetic energy in 
the geographic reference frame), we shall see that this feedback can be positive, 
destabilizing the wobble-less equilibrium. 


Fig. 2 Free precession of a 
deformable Earth. R is the 
instantaneous North Pole 
and 5 is the symmetry pole. 
The geographic coordinates 
(x, y) are centered at O, which 
corresponds to the average 
position of R. 
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3 Dissipation and maintenance 

In the body frame (see Fig. [Til, we write 




( 8 ) 


with = y The magnitude squared of the angular momentum M is 

therefore 

. (9) 

Internal forces keep M fixed, as they produce no net torque. The kinetic energy of 
rotation is 

£ = i(/f2i+/3X2|)=£o + ^(^l-0f2i, (10) 

where Eg = 121^ is the minimum energy consistent with conservation of M. 

Let a = arcsin(f2x/f2) be the angle between Q. and jC 3 , in radians. For small a 
we may express the wobble’s energy as 



For the Earth this gives ^ a} x 10^^ J. 

Displacement of the equatorial bulge (and therefore of the symmetry pole S in 
Fig. 121) induces a small variation in ocean levels, with the same period as the wob¬ 
ble’s. Friction as this “pole tide” moves through straits and shallow oceans dissi¬ 
pates Fw Another source of dissipation is the anelasticity of molten material under 
the Earth’s crust. We would therefore expect the amplitude of the free wobble to 
decay as 

a(f) = (12) 

with Oo set by the Earth’s state of rotation at the time f = 0 when it solidified, and T 
a time constant small compared to the age of the Earth. In terms of the quality factor 


gw 


TCOch 

2 


the power that the wobble dissipates is 


Fw 


E\u Wch 
2 w 


— X 10^0 W . 

gw 


(13) 


(14) 


The fact that a, averaged over many periods of the wobble, remains of order 10^® 
indicates that some mechanism injects into it an average power of order g^^ x 10^ 
W, compensating for the dissipated Pw.lHlEll 

Stochastic perturbations to the Earth’s mass distribution, such as might result 
from seismic events, could cause random displacements of ji: 3 , preventing it from 





6 


Alejandio Jenkins 


becoming aligned with £l. Jeffreys developed a method for estimating gw from the 
data, assuming that the wobble is re-excited stochastically 0. Jeffreys and others 
have thus obtained estimates for gw ranging from 37 to 1,000 Q. These estimates 
are not only uncertain, they generally require more dissipation than can be accounted 
for by known mechanisms of tidal friction and mantle anelasticities 19] [21] |22l, 
while the corresponding Pw seems too high to be compensated by seismic activity 
1^1241 . Moreover, there is statistical evidence that the Chandler wobble is predom¬ 
inantly a deterministic process. 12511^ l27l 

Some authors have therefore concluded that the wobble is forced by geophysical 
fluid circulations ifTTl fl^ [13] [141 . In that case the estimates of gw obtained by Jef- 
freys’s method do not apply directly (a point stressed in ET], but ignored by some 
authors). The key difficulty with such models is that a damped and forced linear 
oscillator ends up moving at the forcing frequency. Even in parametric or non-linear 
resonances, the frequency of steady oscillation is usually a rational multiple ajh of 
the forcing frequency, with integer a and b of low order El. In fact, the polar motion 
has a significant component with a 12-month period, presumably forced by meteo¬ 
rological and other mass transfers connected with the seasons. This seasonal effect 
is subtracted from the data to isolate the Chandler wobble |l3]|2T]|22l. The polar mo¬ 
tion’s power spectrum is inconsistent with a high-quality oscillation (g > 300) being 
driven by an external power spectrum that is smooth in the vicinity of Och ED - The 
observed wobble must therefore be forced by an external power spectrum peaked 
around the resonant Och, which seems like an implausible coincidence. EDEa 
Moreover, it is well known that in the 1920s the amplitude of the Chandler wob¬ 
ble decreased sharply, only to re-start with a phase jump of nearly 180° ED- A 
recent analysis finds two other brief extinctions, followed by phase jumps, in the 
1850s and 2000s IEQI. These irregularities are difficult to explain in stochastic re¬ 
excitation or forced resonance models, as they are not associated with obvious geo¬ 
physical events. 


4 Self-oscillation and intermittence 


Consider the portion of the Earth’s mass made of circulating fluids. In the body 
frame of the solid Earth, these circulations carry significant kinetic energy but neg¬ 
ligible net angular momentum. Resistance to displacing geographically the mean 
route of a circulation may be interpreted as an effective rigidity. With respect to the 
Earth’s rotation, we therefore express the total tensor of inertia as 

= (15) 

with lij given by the distribution of the mass at rest with respect Earth’s solid sur¬ 
face, and lij by flows whose distribution is expressed in Eulerian coordinates (fixed 
geographically), rather than in Lagrangian coordinates associated to individual mass 
elements (see 12 ^ ). 
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The moments of /,y and /,y in Eq. (fTST i are taken to be of the same order, with e 
indicating the small fraction that the circulations contribute to the moments of /J°‘. 
Euler’s equations for the free top may be expressed as 

—+r 2 xM = 0 , (16) 

at 

where d'/dt measures velocity in the body frame JJl. Eor M, = Y.j{hj + 

Eq. (IThl l implies that the equation of motion for the free precession of £l is the 
weighted superposition of the equations of motion for the precessions separately 
induced by Uj and e/jy. 

Let S be the symmetry pole for Ijj. Let us assume, for simplicity, that Uj has its 
own symmetry pole S, as shown in Eig.[3 Eor an asymmetric Ijj, we would consider 
precession about the principal axis of greatest moment of inertia. This precession 
would not be uniform 121 , but our conclusions would be essentially unaffected. 

The centrifugal deformation of Uj propagates with the speed of sound in the 
relevant material, which is very fast compared to the wobble, so that the adjustment 
of S to the displacement of R is nearly instantaneous El. On the other hand, the 
adjustment of S can show an appreciable delay, caused by the circulating fluids’ 
inertia. Eorf? = (x,y), we therefore take 

S{t)=^(x{t-c),y(t-c)) , (17) 

where j3 = OS/OR (Eig.[3]is drawn for constant j3). The precession of R around S 
would then be described by 

ri(f) =-®Eu[y(f)-^y(f-c)] 

\y{t) = d,Enm-l5x{t-c)] ’ 

with ©Eu = ^3 {h ~ I) /I given by the moments of inertia I and for the tensor Itj, 
and by the component O 3 of the angular velocity along the axis of S. In terms of x, 
Eq. (ITSl l corresponds to 

x(f) = —[x(f)+2j3x(f — c) — j3^x(f — 2c)] , (19) 

and equivalently for y. 

Superposing the precession about S, given by Eq. ( fT9] l and weighed by e <C 1, 
on the precession about S given by Eq. El>> taking c <C Tt/cOch (which lets us 
approximate x(f — c) x(f) — cx(f) in the full equation of motion) we obtain 

x + 2cej3(l - j3)G)^yX+®chX = 0 . ( 20 ) 

The coefficient of x in Eq. (l20l i is negative for j3 > 1, indicating that the wobble 
absorbs energy. As long as 


7 = Icep (jS - 1)> coch/Gw 


(21) 
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this effect will dominate over viscous damping, causing R to spiral away from O 
until non-linear dissipative effects match the power input, stabilizing the wobble’s 
amplitude. In the theory of ordinary differential equations, this corresponds to ap¬ 
proaching a weakly non-linear limit cycle; see m . (Note that if the mass associated 
with S did not have significant kinetic energy in the geographic frame, c would be 
negligible and it would be impossible to maintain j3 > 1.) 

The parameters e, j3 vary with the magnitude and distribution of the geophysical 
circulations. The anelasticity of Itj implies that centrifugal deformation would tend 
to make j5 ^ 1, killing the anti-damping in Eq. (l20l i. However, natural changes in 
circulation patterns may counteract this, maintaining j3 > 1 over long periods (recall 
that OR ~ 6 m). 

Instead of tending to a limit cycle, the wobble may start to decay if § approaches 
close enough to the circle OR for the condition of Eq. (l2Tl i to fail. Sensitivity to j3, 
determined by the precise geographical distribution of the circulations, implies that 
extinctions of the wobble need not be associated with major geophysical irregulari¬ 
ties. As OR 0, ^ 1 becomes more likely, and we therefore expect the wobble 

to turn on again, with a random phase jump related to its brief quiescent period. 
The turning on or off of a self-oscillation as the net linear damping changes sign is 
known mathematically as “Hopf bifurcation” (see CH). 

The conditions that the delay c in Eq. (fTTl) remain fixed and very small compared 
to the period of the wobble may be relaxed without qualitatively altering our con¬ 
clusions. Physically, the point is that the wobble itself induces a modulation of the 
geophysical circulations, corresponding to the precession of S in Eig. |3] Eor j3 > 1 
and 0 < OchC < arccos(l/j3), the resulting force on the solid Earth, described by 
Eq. (O, leads the wobble, thus supplying to it net power (see EOl). 


Fig. 3 Torqueless precession 
of a deformable Earth con¬ 
taining fluid circulations. R 
is the instantaneous North 
Pole, 5 is the symmetry pole 
of the portion of the mass at 
rest with respect to the Earth’s 
solid surface, and S is the 
symmetry pole of the part 
of the mass in geophysical 
circulations. 
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5 Outlook 

In the self-oscillatory model offered here, the Chandler wobble is maintained by a 
positive feedback between the wobble and the centrifugal deformation of the part 
of the Earth’s mass in fluid circulations. The energy to maintain the wobble comes 
from the flows, which take it in turn from the heat of solar radiation and of the 
Earth’s internal radioactivity. This is a picture of the wobble, not as an unphysical 
perpetual motion, but as a low-efficiency heat engine, with the geophysical fluids 
playing the part of the working substance and the solid Earth acting as the piston. 
That heat engines and motors in general may be conceptualized as self-oscillators 
has been stressed in ifTSlfT^ . 

Evidence of the modulation of the circulations may already have been detected 
in the atmosphere muiniiii, requiring only re-interpreting that signal as resulting 
from the wobble’s feedback, rather than being an intrinsic property of the circu¬ 
lations. The positive local Lyapunov exponent that ll25l l26l IZTl And in the Chan¬ 
dler wobble data may reflect, not chaos (an aperiodic and strongly non-linear phe¬ 
nomenon), but rather the linear anti-damping of Eq. (|20] | (which makes the wobble¬ 
less equilibrium unstable). The reported variability of this local Lyapunov exponent, 
and its tendency to increase when the El Nino Southern Oscillation (ENSO) is on 
its maximum (El Nino) or minimum (La Nina) temperature phase, seem consistent 
with our model, as ENSO extrema may be associated with stronger atmospheric 
currents around the equatorial region, and therefore with an increased e in Eq. ( fTSl l. 

This simple model raises issues that may be relevant beyond geophysics. A de¬ 
terministic self-oscillation turned on and off by a stochastic parameter may lead 
to behaviors difficult to capture with simple statistical tools. (Eor a review of how 
heavy-tailed distributions may emerge when a stochastic parameter triggers repeated 
dynamical bifurcations, see ED .) This model of the Chandler wobble may therefore 
offer a case study on the importance of better understanding the relations between 
stochastic and deterministic dynamics in complex systems. 
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